GLOBAL POINTWISE ESTIMATES FOR GREEN'S MATRIX OF SECOND 
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Abstract. We establish global pointwise bounds for the Green's matrix for divergence 
form, second order elliptic systems in a domain under the assumption that weak solutions 
of the system vanishing on a portion of the boundary satisfy a certain local boundedness 
estimate. Moreover, we prove that such a local boundedness estimate for weak solutions of 
the system is equivalent to the usual global pointwise bound for the Green's matrix. In the 
scalar case, such an estimate is a consequence of De Giorgi-Moser-Nash theory and holds 
for equations with bounded measurable coefficients in arbitrary domains. In the vectorial 
case, one need to impose certain assumptions on the coefficients of the system as well as 
on domains to obtain such an estimate. We present a unified approach valid for both the 
scalar and vectorial cases and discuss several applications of our result. 



1. Introduction 

In this article, we are concerned with the Green's matrix for elliptic systems 

m m n 

(1.1) Yj L U ui : = ~ Z Z D «( A t D ^)' » = L • • m > 

j=\ j=l a,p=\ 

in a (possibly unbounded) domain Q c R", n > 3. We assume that the coefficients are mea- 
surable functions defined in the whole space R" satisfying the strong ellipticity condition 

m n m n 

(1.2) z Z A y ( ^M^ v ZZi^r =:v ^r> v. e r, 

i,j=l aJ3=l i=l a=l 

and also the uniform boundedness condition 

m n 

(1.3) Yj ZKM ^ Vxer ' 

for some constant v 6 (0, 1]. We do not assume the the coefficients are symmetric. We will 
later impose some further assumptions on the operator but not explicitly on its coefficients. 

In the scalar case (i.e. m = 1), the Green's matrix becomes a scalar function and is 
usually called the Green's function. It is well known that the Green's function G(x, y) is 
nonnegative in D. and for all x, y e Q. with x + y, we have 

(1.4) G(x,y) <K\x-y\ 2 -", 

where K is a constant depending on the dimension n and the ellipticity constant v of the 
operator; see lfl4l . lfT2l . It is also known that if Q. is a bounded domain satisfying the 
uniform exterior cone condition, then for all x, y e D. with x + y, we have 

(1.5) G(x,y) < Kc%\x - y\ 2 ^ a ; d, = disuj, 3D), 
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where K and a e (0, 1) are constants depending only on n, v, and £2; see iTPHl . The methods 
used in iTBll and lfl2ll rely heavily on the Harnack's inequality and the maximum principle 
and does not work for the general vectorial case. By assuming that Q is a bounded C 1 
domain and the coefficients of the operator are uniformly continuous in Q (or belong to 
VMO), Dolzmann and Miiller |6| proved the global estimate ( 11.4b in the vectorial setting. 
It should be noted that Fuchs [7 , 8] obtained a similar result earlier, but under a stronger 
assumption that the coefficients are Holder continuous. Recently, Hofmann and Kim lfl3l 
derived the existence of a Green's matrix in an arbitrary domain, under the assumption 
that weak solutions of the system satisfy interior Holder continuity estimates. They also 
derived various estimates for the Green's matrix of such a system, including an interior 
version of the estimate ( 11.4b . which was applied to the development of the layer potential 
method for equations with complex coefficients in [ 1 1. Their method is interesting because 
it works for both scalar and vectorial cases, but however, they did not attempt to derive the 
global estimates ( 11.4b or ( 11.5b for the Green's matrix in the vectorial setting. 

The goal of this article is to present how one can derive a global estimate correspond- 
ing to (11.4b for Green's matrix of the elliptic systems ( 11.1b in a domain Q. using a local 
boundedness estimate for the weak solutions of the system vanishing on a portion of the 
boundary; see Condition (ILBb below for the precise statement of the local boundedness 
estimate. In fact, we show that such a local boundedness estimate is a necessary and suffi- 
cient condition for the Green's matrix of the system to have a global pointwise bound like 
( 11.41 . We will also show how to derive a global estimate like ( 11.5b for Green's matrix of 
the elliptic system (11.11 ) in a domain Q by using a local Holder continuity estimate for the 
solutions of the system vanishing on a portion of the boundary <9Q; see Condition ( ILHb 
below for the statement of the local Holder estimate, and also the condition (iLH'b in Re- 
mark B.l II which is a little bit weaker. The novelty of our work is in presenting a unifying 
method that re-proves the global estimates ( 11.4b and ( 11.51 ) for the Green's function for 
the uniformly elliptic operators with bounded measurable coefficients as well as the cor- 
responding estimates for the Green's matrix of the elliptic systems (11.1b . for instance, in 
a bounded C 1 domain with uniformly continuous or VMO coefficients. Moreover, it has 
other interesting applications to L°°-perturbation of diagonal systems in a domain satisfy- 
ing the uniform exterior cone condition, elliptic systems satisfying Legendre-Hadamard 
condition in a bounded C 1 domain with principal coefficients in VMO and lower order 
terms in L°°, Stokes systems in a three-dimensional Lipschitz domain, etc.; see Section [4] 
below. As a matter of fact, application to L^-perturbation of diagonal systems in a do- 
main satisfying the uniform exterior cone condition shows the power and flexibility of our 
method since that result does not seem to follow from other known methods, such as that 
based on the LP theory by Dolzmann and Miiller |6|. 

The organization of the paper is as follows. In SectionffJ we introduce some notations 
and definitions including our definition of the Green's matrix of the system (11. lb in Q. 
In Section [3] we give precise statement of the conditions concerning the local estimates 
for weak solutions of the systems and state our main theorems. In Section HI we present 
applications of our main results. The proofs of our main results are given in Section|5]and 
a technical lemma is proved in Appendix. 

Finally, a few remarks are in order. We do not treat the case n — 2 in our paper. In 
two dimension, the Green's matrix has logarithmic growth and requires some different 
methods; see J6J and also 0. As a matter of fact, the method used in this paper breaks 
down and does not work in that case. By this reason, the two dimensional case will be 
discussed in a separate paper (4), where we will also treat a parabolic extension of our 



GREEN'S MATRIX 



3 



results. As alluded earlier, the main difference between our result and that of [13] is that 
they were mostly concerned with the Green's matrix of a L^-perturbed diagonal systems 
in the whole space R" and focused on interior estimates of the Green's matrix while our 
paper is mainly concerned with the global estimates like ( 11.41 ) and ( 11.5b . which we believe 
are quite more useful in practice, especially in the vectorial case. In 0, Auscher and 
Tchamitchian introduced the "Dirichlet property (D)" in connection with the Gaussian 
estimates for the heat kernel of the operator L, which is very similar to the condition (ILHI l 
of this article. We would like to hereby thank Pascal Auscher for kindly informing us about 
the paper 

2. Notations and Definitions 

Let L be an elliptic operator acting on column vector valued functions u — (u l , . . . , u"') T 
defined on a domain Q c R", n > 3, in the following way: 

Lu = -D a {A ap D p u), 

where we use the usual summation convention over repeated indices a,B = 1, . . . ,n, and 
A a P = A a/1 (x) are m x m matrix valued functions on R" with entries A"? that satisfy the con- 
ditions ( 11.2b and d 1 . 3b - Notice that the ;-th component of the column vector Lu coincides 
with Yij Lijitj in ( 11.1b . The adjoint operator L of L is defined by 

'Lu = -D^A^Dpu), 

where A a ? = (A^ a ) T ; i.e., 'A"f = Af. We use the same function space Y h2 (Cl) as in {HI . 
For reader's convenience, we reproduce the definition below. 

Definition 2.1. For an open set Q. c R" (« > 3), the space y 12 (Q) is defined as the 
family of all weakly differentiable functions u e L 2 "/ ( "~ 2) (Q), whose weak derivatives are 
functions in L 2 (£2). The space Y l,2 (Q.) is endowed with the norm 

IMIy«(n) := IMb'/Mfn) + I|Dm|| L 2 (£2) . 

We define Fq' 2 (Q) as the closure of C™(Q.) in F U (Q), where C™(Q) is the set of all infin- 
itely differentiable functions with compact supports in fl 

Remark 2.2. If |Q| < oo, then Holder's inequality implies F 1,2 (f2) c W l,2 (£l). In the case 
Q. = R", we have F U (R") = ^(R"). Notice that by the Sobolev inequality, it follows that 

(2.3) Nlzwp) < C(n)||D M || L 2 (n)) V« 6 Y l Q ' 2 (Q). 

Therefore, we have W^iSl) c Y^iSl) and W^ 2 {Q.) = K ( J ,2 (Q) if |0| < oo; see fT5l . 

Definition 2.4. Let 2" c Q. and u be a Y l,2 (Q) function. We say that u vanishes (or write 
u — 0) on E if u is a limit in F 12 (Q) of a sequence of functions in C t M (Q \ £). 

Notation 2.5. We denote by L™(Q) the family of all L°°(Q) functions with compact sup- 
ports in £1. Notice that L™(Q) = L°°(Q) if Q is bounded. 

Notation 2.6. We denote D. R (x) = D. n B R (x) and Z R (x) - dD. n B R (x) for any R > 0. We 
abbreviate Q R = Cl R {x) and £ R = E R {x) if the point x is well understood in the context. 

Definition 2.7. We say that an m x m matrix valued function G(x,y), with entries Gij(x,y) 
defined on the set {(x,y) e Q x £2 : x + y}, is a Green's matrix of L in Q if it satisfies the 
following properties: 
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i) G(-,y) € WMCQ) and LG(-,y) = 6 y I for ally 6 Q, in the sense that 



(2.8) f A a ^DpGjk{-,y)D a (p i = <p k (y), = (0 1 . . . ,cf> m ) T e C C °°(Q). 



J<2 



ii) G(-,y) £ T U (Q \ B r (y)) for all y e Q and r > and G(-,y) vanishes on <9Q. 

iii) For any / = (/',.. .,f m ) T e L™(Q), the function u given by 




belongs to Yb 2 (Q) and satisfies 'Lu — f in the sense that 



(2.10) f A°jD a u'Dp<y = f /V, V0 = (0 1 ,...,0" , ) r eC c oo (O). 



Jn Jn 



We note that part iii) of the above definition gives the uniqueness of a Green's matrix; 
see (T3). We shall hereafter say that G(x,y) is "the" Green's matrix of L in D. if it satisfies 
all the above properties. 



The following condition, which hereafter shall be referred to as dLBt . is used to obtain 
pointwise bounds for the Green's matrix G(x,y) of L in CI. 

Condition (LB). There exist R max e (0, oo] and No > such that for all x e £1, R e 
(0, R max ), and / e L°°(Q s (x)), the following holds: If u e W u (Q R (x)) is a weak solution 
of either Lu - / or 'Lu — f in Q«(x) and vanishes on £r(x), then we have 



Notation 3.1. We use the convention c ■ oo = oo for c > and l/oo = 0. 

Remark 3.2. By using a standard covering argument, it is easy to see that the constant Rmax 
in the condition ( ILBI i is interchangeable with c-R max for any fixed c e (0, oo), possibly at the 
cost of changing the constant No in the condition (ILBb by • A^o, where K = /f(n, c) > 0. 

Theorem 3.3. Assume the condition ( ILBt a«c/ /ef G(x,y) be the Green's matrix of L in CI. 
Then we have 

(3.4) \G(x, y)\ < C\x - yt" for all x,ye£l satisfying < \x - y\ < R„ MX , 

where C — C(n, m, v, iVo). 

The following condition (IIHI i combines two conditions that appeared as the properties 
(H) and (H); oc in JT3). It means that weak solutions of Lu =0 and 'Lu — in B c Q. are 
locally Holder continuous in B with an exponent jjo- 

Notation 3.5. We denote aAb = min(a, b) and d x = dist(x, dQ). 

Condition (IH). There exist // e (0,1], R c e (0, oo], and N\ > such that for all x e CI 
and R € (0,d x A R c ), the following holds: If u € W ia (B R (x)) is a weak solution of either 
Lu = or 'Lu = in Br(x), then we have 



3. Main results 



(LB) 



Hh|Il~(£V) 



< No (R-" /2 



?(fi R) +R 2 \\f\\L-(n R )); CIr:=C1r(x). 



(IH) 




for < r < s < R. 
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Theorem 3.6. Assume the conditions (lIHb and (ILBt . Then, the Green 's matrix G(x, y) ofL 
in Q exists and satisfies the estimate ( 13.41 ) with C — C(n, m, v, No). Also, the Green 's matrix 
'G(x,y) of the adjoint operator L in Q exists and we have 

(3.7) G(x,y) = 'G(y,x) T , V.tjed, x + y. 
Moreover, the Green 's matrix G(x, y) satisfies the estimate 

(3.8) \\G(;y)\\ Pm \B r( y)) < C/ 2 -" )/2 , Vre(0,/W), Vy € £2, 
where C — C(n, m, v, No). 

The following theorem says that the converse of Theorem |3.3| is also true if we assume 
the condition (lIHI i. 

Theorem 3.9. Assume the condition (lIHI i and let G(x,y) be the Green's matrix ofL in Q. 
Suppose there exists R max € (0, oo] such that for all x, y € Q. satisfying < \x — y\ < R max , 
we have 

(3.10) \G(x,y)\<Co\x-y\ 2 -". 

Then the condition dLBl l is satisfied with the same R max and No — N)(n, m, v, Co). 

Remark 3.11. In fact, one can replace the condition (lIHI i in Theorem l3.6l and Theorem |3.9l 
by the following condition ( lIH't : There exist po e (0, 1], R c e (0, oo], and Co > such that 
if u e W x ' 2 {B R {x)) is a weak solution of either Lu — or 'Lu — in Br(x), where x e £2 
and < R < d x A R c , then we have 

(IH') ["W^w^Co/r^f |h| 2 

Vf! fl (.v) 

Here, [n]po(%) denotes the usual Holder seminorm. It is not hard to see that condition 
(lIHI i implies (lIH'b with Co = Co(n, m, v, //o, M ) and the same po and R c . As a matter of fact, 
the conditions dlHI ) and dlHl are equivalent under our basic assumptions on the operator 
L; see |[T3l Lemma 2.3]. However, the condition (IIHI does not imply dlHt , for instance, 
in the presence of lower order terms in the operator L. In this sense, UH'i is a weaker 
condition. We point out that the properties (H) and (H); oc in lfT3l can be replaced entirely 
by the condition dIH'l) . without affecting the conclusion of the main theorems in lfT3l . 

The following condition, which hereafter shall be referred to as dLHI ). is a combination 
of dIHb and another condition that appeared as the property (BH) in lfl3l . 

Condition (LH). There exist fx € (0, 1], R„ mx e (0, oo], and N\ > such that for all x € Q 
and R e (0,R max ), the following holds: If u 6 W ia (Q. R (x)) is a weak solution of either 
Lu = or 'Lu = in Q.r(x) and vanishes on £r(x), then we have 

(LH) f \Du\ 2 < N\ (-)" 2+2/J ° f |Z)h| 2 for < r < s < R. 

Remark 3.12. In the condition (ILHb . the constant R max is interchangeable with c • 7? mflJC for 
any fixed c e (0, oo), possibly at the cost of changing the constant N\ in the condition (ILHI l 
by K ■ Nu where = K(n, c) > 0. 

It will be shown in Appendix that dLHb implies dLBI i. Also, it is obvious that (ILHI l 
implies dTFTb. Therefore if (ILHb is satisfied, then by Theorem 13.61 the Green's matrix 
G(x,y) of L in Q exists and satisfies the estimate (13.4b . The following theorem asserts that 
in fact, in such a case, a better estimate for G(x, y) is available near the boundary dQ.. 
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Theorem 3.13. Assume the condition (ILHI l and let G(x,y) be the Green's matrix of L in 
Q,. Then for all ijefi satisfying < \x — y\ < Rmax, we have 

(3.14) \G(x,y)\ < C\d x A \x-y\r[d y A \ x - y \f>\ x - y f-^\ 

where C — C(n,m,v,/j.o,Ni). If Rmax < 00 an d £2 is bounded, then for all x,y 6 Q with 
x + y, we have the estimate ( 13.14b with C — C(n, m, v, //o, Ni , Rmax/ diam £2). 

Remark 3.15. It will be clear from the proof of Theorem l3 . 1 31 (see ( 15.391 ) in 35.41 ) that we 
in fact have the following estimate, which is slightly stronger than ( 13.141 ): For all x,y e Q 
with x + y, we have 

\G(x,y)\ < C\d x \\x-y\A R max r\d y /\\x — y\/\ RmaxT{\x - yl A R^'^r 

where C = C(n,m,v,fio,N\). 

Remark 3.16. The following condition ( ILHI can be used as a substitute for ( ILHI ) in The- 
orem [3T3] There exist yUo e (0,1], Rmax 6 (0, °°], and Co > such that for all x e CI and 
6 (0,Rmax), the following holds: If u e W x ' 2 {Q. R (x)) is a weak solution of either Lu — 
or 'Lh = in Q«(x) and vanishes on E R {x), then we have 

(LH') mcm(B m {x)) <C R-^(-f \u\ 2 ) , where u =*n flW «. 

It is not hard to verify that the condition ( ILHI l implies ( ILHI with Co = Co(«, m, v, yuo, iV] ) 
and the same fig and Rmax- Also, it can be easily seen that the condition ( ILHI implies both 
the conditions dLBI ) and dlH't . We note that the condition dLH'b is, however, weaker than 
condition (ILHI l in general. From the proof of Theorem 13. 131 it should be clear that the 
conclusion of Theorem l3. 13l remains the same under a weaker condition (ILH'b . 

4. Applications of Main Results 

4.1. Scalar case. In the scalar case (i.e., m = 1), both conditions ( ILBt and (lIHt are sat- 
isfied with Rmax - 00 an d A^o = No(n,v) in any domain O; see e.g., IfTTI . Also, in the 
scalar case, the Green's matrices are nonnegative scalar functions; see lfT2l . Therefore, the 
following corollary is an immediate consequence of Theorem l3.6l 

Corollary 4.1. Ifm — 1, then for any domain Q, c R", the Green 's function G(x,y) of L in 
Q. exists and satisfies 

(4.2) G(x,y) < C\x-y\ 2 -'\ Vjc,yeQ, x + y, 

where C — C(n, v) is a universal constant independent of £1. 

Remark 4.3. Corollary 14. 1 1 is widely known (see e.g., Ifl2l \14\ ). However, it should be 
mentioned that unlike fl2l [T4"l . we do not need to assume that Q is bounded. 

Also, in the scalar case, the condition ( ILHI ) is satisfied if Q satisfies the condition ([S]), 
the definition of which is given below. In fact, if L is a small L°°-perturbation of a diagonal 
system, then the condition ( ILHI ) is satisfied whenever Q satisfies the condition <JSj ; see 34.21 
below. 

Condition (S). There exist 6 > and R a e (0, 00] such that 

(S) \B R (x)\a\>6\B R (x)\, Vxe«9Q, \fR&(0,R a ). 



The following corollary is then an easy consequence of Theorem |3.13l 
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Corollary 4.4. Assume m — 1 and let G(x, y) be the Green 's function ofL in Q., where £2 is 
a domain satisfying the condition |[S). Then, G(x, y) satisfies the estimate ( 14.2b . Moreover, 
for all x, y € Q satisfying < \x — y\ < R a , we have 

(4.5) G(x,y) < C[d x A \x-y\f°{d y A \x - y\f°\x - y \ 2 -"- 2 '">, 

where C — C(n, v, 6) and fiQ — /Uo(", v, 6). IfR a < 00 and Q. is bounded, then for all x, y e Q 
with x +y,we have the estimate (14. 5t with C — Cin, v, 0, R a / diam Q). 

Example 4.6. £1 = R" satisfies the conditions (© with 0=1/2 and R a = oo. Therefore, 
Corollary |4.4| implies that for all x,y € R" with x + y, we have 

G(x,y) < C{d x A |x — yin^ v A \x - y\f°\x - y \ 2 -"- 2 ''\ 

where C = C(n, v) and //o = jUo(«, v). 

4.2. L°° -perturbation of diagonal systems. Let a a P(x) be scalar functions satisfying 

n 

(4.7) a a P{x)^ a >v*\tf, V^eR"; £ |a^(x)| 2 < Vq 2 , 

for all x E R" with some constant vo € (0, 1]. Assume that Q satisfies the condition ([S} and 
let A"f(x) be the coefficients of the operator L. We denote 



(4.8) 




where dij is the usual Kronecker delta symbol. By j 13, Lemma 4.6], there exists a number 
Sq = <§o(n, vq, 6) such that if S < Sq, then the condition (ILHI) is satisfied by L in £1 with 
parameters fj.Q = /uo{n, vq, 0), N\ = N\{n, m, vo, 6), and R max = R a - Therefore, the following 
corollary is another easy consequence of Theorem l3.6l and Theorem l3. 131 

Corollary 4.9. Let a a @(x) satisfy the condition ( 14.7b . Assume that Q. satisfies the condition 
<[S} an d let S be defined as in ( 14.81 ), where A"j(x) are the coefficients of the operator L. 
There exists Sq — So(n, vq, 6) such that if § < Sq, then the Green 's matrix G(x, y) ofL in Q 
exists and for all x, y e Q. satisfying < |x — y\ < R a , we have 

(4.10) \G(x,y)\ < C{d x A \x - y\f a {d y A \x - y\f°\x - y \ 2 -"- 2 ^, 

where C — C(n,m,vo,6) and fXo — //()(«, Vo, &)■ If R a < 00 and Q is bounded, then for all 
x, y 6 Q such that x +y,we have the estimate ( 14.10b with C — C(n, m, vo, 0, R c ,l diam Q). 

Example 4.11. Let fi = jie R" : x„ > if(x')}, where x = (x', x„) and <p : R" _1 -> R is a 
Lipschitz function with the Lipschitz constant K. Then Q satisfies the condition (O with 
6 — 0(n, K) and R a — oo. If L is a small L°° -perturbation of a diagonal system in the sense 
of Corollary 14. 91 then the Green's matrix G(x,y) of L in Q exists and we have 

\G(x,y)\ < C{d x A |x - y\f°{d y A \x - y\f a \x - y| 2 ~"~ 2/ \ Vx, y eQ., x+ y, 

where C = C(n, m, vq, K) and /jq = p.o(n, vq, K). 
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4.3. Systems with VMO coefficients. For a measurable function / on R", we set 

to s (f) := sup sup + \f(y) - f x , r \ dy, Vc5 > 0; f„ =\ f. 

xeW r<5 J B,.(x) J B r (x) 

We say that / belongs to VMO if lim^o oj s {f) = 0; see 02). 

If the coefficients A a P of the operator L are functions in VMO satisfying dl.2t . dl.3| > 
and if O is a bounded C 1 domain, then the condition dLHI ) is satisfied with parameters po, 
Ni, and R max depending on O and a>g(A a ^) as well as on n,m,v. Therefore, we have the 
following corollary of Theorem l3.13l 

Corollary 4.12. Let Q be a bounded C l domain. Suppose the coefficients A a @ of the 
operator L belong to VMO and satisfy the conditions ( 11.2b , (11.3b . Then for all x, y E Q 
with x + y,we have 

\G(x,y)\ < C\d x A |jc - y\r{d y A\x- y|p|x - yf-"^", 

where C and po are constants depending on n, m, v, Q, and a) 6 (A a/3 ). 

In the above corollary, one may assume that A a/S satisfy the weaker Legendre-Hadamard 
condition and may even include lower order terms in the operator. More precisely, let 

(4.13) L A u = -D a {A al3 D B u) + D a (B a u) + B a D a u + Cu + Au, 

where A a/3 , B a , B a , and C are m x m matrix valued functions on R" satisfying 

A"f(x)^^%r] a > v|f| 2 |i/| 2 , V£ E R ffl , V?7 6 R", Vx E R"; 

(4 - 14) zirt^v- z(irt + rt) + iict,v- 

a,B=l a=\ 

for some constant v E (0, 1], and A is a scalar constant. 

Corollary 4.15. Let Q be a bounded C 1 domain and let the operator Lx be defined as 
in d4.13t with the coefficients satisfying the conditions ( 14. 14b . We assume further that the 
leading coefficients A a P belong to VMO. There exists Aq > such that if A > Aq, then the 
Green 's matrix G(x, y) ofLx in Q exists and for all x, y E Q with x +y,we have 

\G(x,y)\ < C{d x A \x - y\r\d y A \x - y\T\x - yf-"- 2 ^, 

where the constants po and C depend on n, m, v, Q, A, and u>s(A a P). 

To give a sketch of proof for Corollary 14. 151 first we note that for sufficiently large A, 

12 12 

we have the solvability of the following problem in Y () ' (Q) m = W '~(Q)"': 

J L,\u - f in Q, 
[ u — on d£l, 

where / E L)f (Q). In particular, one can construct the "averaged Green's matrix" G p (x,y) 
of Lx in £2 by following the argument in Ifl3l §4]. Also, it is not hard to see that the 
condition dLHI in Remark D.l H is satisfied in this case. In particular, we have the condition 
dLBb . We point out that these are all the ingredients needed for construction of the Green's 
matrix G(x, y) of Lx in £1 Then by modifying the proofs of Theorem |3.6| and Theorem |3.131 
one can prove the above corollary; see Remark [3.16l The details are left the the reader. 

Remark 4.16. In Corollary 14. 121 and Corollary 14. 151 the conditions of Q and A"* 3 can be 
relaxed. We may assume that Q is a bounded Lipschitz domain with a sufficiently small 
Lipschitz constant, and a>s(A a P) is also sufficiently small for some 6 > 0; see e.g., J2). 
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4.4. Stationary Stokes system. Let Q c R 3 be a bounded Lipschitz domain with con- 
nected boundary. We consider the stationary Stokes system 

(4.17) -Ah + V/? = 0, div«=0 in Q. 

It is known that the condition (ILHI l is satisfied in this setting; see ifTTll and also 0. We 
also note that Caccioppoli's inequalities are available for the system ( 14. 17b . Then again, 
by modifying the proof of Theorem l3.13l one can prove the following corollary. 

Corollary 4.18. Let Q c R 3 be a bounded Lipschitz domain with connected boundary. Let 
G be the Green's matrix of the stationary Stokes system J4. 17b in il Then for all x,y e Q, 
with x +y,we have 

(4.19) \G(x,y)\ < C{d x A \x-y\f{d y A \x - y\f°\x - y\- l ~ 2 ^ 

for some positive constants C and fXo depending on Q. 

We remark that estimate ( 14.191 ) of the Green's matrices for the Lame system and the 
Stokes system were recently shown in J3) by a different method to ours. 

5. Proofs of Main Theorems 

5.1. Proof of Theorem 13.31 Let R e (0, R max ) and y e Q be arbitrary, but fixed. Assume 
that / e L°°(Q) is supported in Qs(y) and let u be defined by ( 12.91 . Notice that we may 
take u in place of in ( 12.101 . Then by ( 12.31 we have 

(5.1) \\u\\ L 2„n„-2 )(n) < C\\Du\\ LHn) < C||/|| L 2»/ ( „ +2 , (n) < C/? 1+ " /2 ||/|| L » ( n fit v))- 

Also, notice from Remark F2~2l that u e W l ' 2 (Q.R(y)). Therefore, u is a weak solution of 
Lu — f in Cluiy) vanishing on 2«(j) and thus, by the condition (ILBI) we have 

(5.2) NMowW) * ^Vo (R-" /2 \\u\\ L 2 (nM) + ^ 2 ||/llL» ( n R tv))) • 
Then by ( 15.21 . ( 15.11 . and Holder's inequality, we derive 

(5.3) IMIz.°°(n s/2 (y)) < C/? 2 ||/||L»(n s ( V )). 
Hence, by (12.9b and ( 15.31 . we conclude that 

(5.4) f G(-,y)f 

Ja R (y) 

Therefore, by duality, we conclude from d5.4t that 

(5.5) IIGC.y)lbm,cy» ^ cr1 - 

Next, notice that (ILBI implies that for x e O and R e (0,R max ), if v e W l - 2 (Q. R (x)) is a 
weak solution of Lv — in Q.r = Q.r{x) vanishing on Zr{x), then we have 

I|v||l» ( q s/2 ) < N R-" p -\\v\\ L 2 ( n R) . 
Then, by a standard argument (see e.g., |[T0l pp. 80-82]) we also have 

(5.6) llv|k~ ( n s/2 ) < C p R-""'\\v\\ L „ (nRh Vp > 0, 

where the constant C p depends on n, No, and p. 

Now, for any x € £2 such that < \x - y\ < Rmax/2, set R := 2\x - y|/3. Notice that 
Definition 12 . 7 l implies that G(-,y) e W 1,2 (Q.r(x)) and satisfies LG(-,y) = weakly in Qr(x) 
and G(-,y) = on Er(x). Therefore, by ( 15.61 and ( 15.51 , we have 

(5.7) \G(x,y)\ < CR-"\\G(;y)\\ V(nR(x)) < CR-"\\G(;y)\\ V{niRiv)) < CR 2 ". 



< CR 2 \\f\\ L ~ { n R{y)h V/ e L°°(Cl R (y)). 
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We have thus shown that 

\G(x, y)\ < C\x - y\ 2 ~" for all x, y e O satisfying < \x - y\ < R ma x/2, 
where C = C(n, m, v, No). The theorem then follows from Remark l3~2l ■ 

5.2. Proof of Theorem l3.6l The existence of the Green's matrices G(x,y) and 'G(x,y) as 
well as the identity (13.7b is a consequence of the condition ( lIHk see ifTJl Theorem 4. 1 ] and 
lfT3l Eq. (4.34)]. Then, by TheoremES the condition dLBl yields the estimate d3~4l . Also, 
by El Eq. (4.24)], we find that the estimate in (O is valid for < r < (d y A R c )/2. To 
give a full proof of ( 13.81 l. we need to make use of the estimate (13.41 ) and adapt the arguments 
used in IfTTI as follows. 

For p > 0, let C^y) e 7o' 2 (Q) be the averaged Green's matrix of L in D. as constructed 
in El §4.1]. Notice that by H3 Eq. (4.3)], we have 

(5.8) f kfD^A-^Dj = f «*, Vh € F^O). 

Jo. Jn p (y) 

Also, by HI Eq. (4.2)], we have 

(5.9) UDCi-^hm) ^ QSl p {y)f- n)l2n < Cp (2 -" )/2 , 

1 2 

where C = C(n, m, v). Denote by H the Hilbert space Y Q ' (Q) m with the inner product 

<m,v) := I D a u l D a v\ 
Jn 

For all / € L~(Q), the linear functional 

C fw 



Hi 

is bounded on H. Hence, by the Lax-Milgram lemma there is a unique u e H satisfying 



(5.10) f AlfD p w j D a ii = f/-H>, VweH. 

Thus, if we set w to be the £-th column of G p {-,y) in ( I5.10l l, we obtain from ( 15.81 ) that 

(5.H) f G£(-,;y)/ < = f «*. 

Jn Jn P (y) 

Also, if we set w — u in ( I5.10l l, it follows from ( 12.3b that 

(5.12) IMIzW ( n) < C||Dit|b (n) < C\\f\\o-i^(ay 

Let us now assume that / is supported in Q.r :— Q«(y), where y e Q. and R e (Q,R max ) 
are arbitrary, but fixed. Notice that u e W 1,2 (f2«) and h is a weak solution of 'Lu = f in Q.r 
vanishing on Sr. Therefore, the condition ( ILBI ) implies that 

NIl-OTw) < #o (*-" /2 IMb ( n s ) + ^ 2 H/llL» ( n s )) • 
On the other hand, ( 15.12b and Holder's inequality yields 

Mb<n.) < C/? 2+ " /2 ||/||z.» ( n s ) 
By combining the above two inequalities, we obtain 

(5.13) ll«llL» ( n„ /2 , < CR 2 \\f\\ L ~ ( n R) . 
Then, by ( 15.111 ) and ( 15.131 ) we derive 



f 

Jn„ 



G^(-.y)/ ; 



< C/? 2 ||/|| L ~ ( n R ), V/ e L°°(C1 R ), Vp e (0, R/2). 
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Therefore, by duality, we conclude that 

W&i-^WLHQxW^CR 2 , Vpe(0,*/2). 

Now, for any x e Q such that < \x - y\ < R max /2, let us take R := 2\x - y\/3. Notice that 
if p < R/2, then G p (-,^) e W h2 (D. R (x)) and satisfies W(-,y) = in fi fl (jc) and vanishes 
on 2r(x). Therefore, as in ( 15.71 , we have 

< C/?-"||G^(-,y)|| Ll( n 3s(v)) < CR 2 ". 

We have thus proved that for any x,y e Q. satisfying < \x - y\ < R max /2, we have 

(5.14) Wix^KQx-yt", Vp<\x-y\/3, 

where C = C(n, m, v, No). 

Next, fix any r e (0, R max /2) and let v p be the k-th column of G p (-, y), where k — l,...,m 
and < p < r/6. Let 77 be a smooth function on R" satisfying 

(5.15) < 77 < 1, r) = 1 on R" \ B r (y), 77 = on B r / 2 00, and ID77I < 4/r. 
We set 11 = j] 2 v p in (15.8b and then use (15.14b to obtain 



(5.16) f 77 2 |£>v p | 2 <C f |Dt7| 2 |v p | 2 < Cr" 2 f |x - yl 2 ^"' da: < Cr 

Jn Jn J B,(y)\B r / 2 (y) 

Therefore, by (15.15b . (12.3b . and ( 15.16b , we obtain 



provided that < p < r/6. On the other hand, if p > r/6, then ( 15.9b implies 

lklU-> ( Q^ W) * IklU^w < c|K|| i2(n) < c/ 2 -«>/ 2 . 

By combining the above two estimates, we obtain 

(5.17) ||G"(-,y)|lL^)(n\B r (y))<C/ 2 -" )/2 , Vre(0,/W/2), Vp > 0. 
Notice from (15.16b and (15.15b that for < p < r/6, we have 

\\DGP(;y)\\ LW(y)) <CrV- n)/2 . 
In the case when p > r/6, we obtain from ( 15.9b that 

l|DC(-,:>0lb(n\ flrW ) < \\DGP(;y)\\ LHn) < Cp^' 2 < Cr {2 - n)l2 . 
By combining the above two inequalities, we obtain 

(5.18) \\DG p {-,y)\\ L i ( n\B, m <Cr {2 - n) l 2 , Vr e (0,R,„ ax /2), Vp > 0. 

Notice from lfl3l Eq. (4.19)] that there exists a sequence {p P j^° =1 tending to zero such 

that CH-.y) ->■ G(-,y) weakly in F Q 1,2 (fi \ B r (y)) for all r > 0. Therefore, <|3JD follows 
from ( 15.17b . ( 15.18b . and the obvious fact that R max /2 and R max are comparable to each other 
in the case when R max < 00. The proof is complete. ■ 
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5.3. Proof of Theorem l3.9l As we mentioned in the proof of Theorem l3.6l the condition 
(lIHb implies the existence of the Green's matrices G(x,y) and 'G(x,y) in Q and also the 
identity (13.7l i. Hence, if G(x,y) satisfies the estimate (13 . 1 Ob . then so does 'G(x,y). There- 
fore, by the symmetry, it is enough to prove (ILBI i for weak solutions of 'Lu = /. 

Let x € Q and R e (0,R max ) be given. Assume that u e W 1 ' 2 (Q.r(x)) is a weak solution 
of 'Lu = / in Q.r(x) vanishing oiiZr(x), where / € L°°(f2s(x)). Then, we have 

(5.19) f A^D a u%w j = f /V, Vh- e W ( ;- 2 (Q«W). 
Jn„0) 

Let G p (-,y) be the averaged Green's matrix of L in Q as in the proof of Theorem 13. 61 
Set v = £u, where £ is a smooth cut-off function on R" satisfying 

(5.20) 0<£<1, supp£ c B R/2 (x), £ = 1 on B 3R/& (x), and < 16/7f. 
Notice that v e rd' 2 (Q) and thus ' bv ED , we obtain 

(5.21) -f = f A^D p G p jk {;y)u i D a C+ f A"jDpG p k (; y)£D a u'. 

Jo p (y) Jo Jo 

On the other hand, notice that £G p (;y) e Wi^n^C*:)). Hence, if we set w to be the A:-th 
column of (G p {-,y), then by ( 15.191 l, we obtain 

(5.22) f A^D a u i G p k (;y)Dtf+ f A^D a u^D fi G p (;y) = f (f j G p (-,y). 

Jo R (x) Jq„oo Jn fl (.v) 

Recall that supp£ c Br^x). Therefore, by combining (15. 2U and (15.22t , we obtain 

(5.23) -f £u k = f A^D p G p jk {;y)u i D a C- f A^&.^yWDtf + f Cf j G p k (-,y) 

jQ. p (y) JO. JQ JO. 

=: /] +/ 2 + / 3 . 

Now, assume that y € Q«/4(jc). Notice from d5.20| i that we have dist(y, supp D() > /?/8. 
Set r — R/8 A (c/ v A R c ). By flT3] Eq. (4.17)] and (BJ Eq. (4.19)], there exists a sequence 
IP/i}~ =1 tending to zero such that G p, '(-,y) ->■ G(-,y) weakly in W hq (B r (y)) for q e (1, 
and G p *(->)0 -» G(-,y) weakly in Y^ 2 (Q. \ B r (y)). Notice that 

h+h= f A :^ DpG p k {-,y)u'D a % - f A"?G p J-,y)D a u i D^\ 

Jo\B r (y) Jo\B r (y) 

h= f tf j G p (;y)+ f {f j G p (;y). 

Jn\B r (y) JB r (y) 

Therefore, by taking limits in ( 15.231 1 and using d5.20| i, we have for almost all y e Q«/4(x), 

(5.24) u k (y)= f AfiDpG^yWDcC- f A°fG jk (;y)D aU i D B { + f {f j G jk (;y) 

Jo Jo Jo 

=:l[+l' 2 +l' v 

Denote A«(y) = Q.3R/i(y) \ Bu/^(y). By using Holder's inequality and d5.20| ) we obtain 

< C^ 1 ||Z)G(.,y)llz?CA^))Nb(B w W). 
\4<CR- l \\G(;y)\\ \\Du\\ 

and 

|4| < C||G(-,y)|| L . (n3S/4(v)) ||/|| L » ( n s/2 (.v)). 



GREEN'S MATRIX 



13 



Denote A R (y) = Cl R (y) \ B R/ i 6 (y). Observe that G(-,y) e W h2 (A R (y)) and it satisfies 
LG{-,y) = weakly in A R (y) and vanishes on d£l D dA R (y). Then by the Caccioppoli's 
inequality and the estimate (13 . 1 Ob . we obtain 

\\DG{;y)\\ LHAR(y)) < C J R- 1 ||G(-,y)|| L2(Xstv)) < CR^" )I2 . 

Therefore, we have 

(5.25) < Ci?-" /2 || M |b (nsW) . 

By setting w = ifu in ( 15.191 l, where r\ e C™(B R (x)) is a cut-off function such that /jslon 
Br/2(x) and ID77I < 4-/R, and use a standard argument, we derive 

f i] 2 \Du\ 2 <C f \D n \ 2 \u\ 2 + C f \T}f\\TlU\. 

Jn R (x) Jn R (x) Jn R (x) 



By the Sobolev inequality, Holder's inequality, and Cauchy's inequality, we obtain 

Jn R (x) 

< e 



f 



f |Dt,| 2 |m| 2 + e f // 2 |Dh| 2 + Cir 1 ^ +2 ||/|£.~ W) . 

By choosing e small enough, we then obtain 

f i] 2 \Du\ 2 <C f \Dtf\u\ 2 + CR" +2 \\f\^ {ClR{x)) . 
Ja R (x) Jn R (x) 

Therefore, by using the estimate ( 13.10b we derive 

(5.26) \l' 2 \ < CR- n/2 \\u\\ LH£lR(x)) + CR 2 \\fh~ { n Ri x)). 
By using the estimate ( 13.10b again, we also obtain 

(5.27) \I' 3 \ < CR 2 \\f\\ L ~ { n R{ x)). 

By combining above estimates ( 15.251 ), ( 15.261 ), and ( 15.27b , we conclude from ( 15.24b that 

(5.28) NMqwM) < C (R-" /2 \\u\\ LHns) + R 2 \\f\\ L ^ R) ) , 

where C = C{n,m,v,Co). Since ( 15.28b holds for all x e Q. and R e (Q,R max ), we obtain 
dLBb by a standard covering argument. ■ 

5.4. Proof of Theorem l3.13l Notice that by Lemma l6.1l and Theorem l3.61 we have 

(5.29) |G(x,.y)|<C„|x-3;| 2 -" if < \x - y\ < R max , 

where Co = Co(«, m, v,jUo, N\). To prove the estimate ( 13.14b . we first claim that 

(5.30) \G(x, y)\ < C{d x A \x - y\f a \x - yt"^ if < \x -y\< R max , 
where C = C(n, m, v, [io, Ni). The following lemma is the key to prove (15.30b . 

Lemma 5.31. Assume the condition (1LHI) . For R e (0, R ma x) and x e £1 such that d x < R/2, 
letu e W l ' 2 (Q. R (x)) be a weak solution of Lu = in £l R (x) vanishing on Z R (x). Then, we 
have 

(5.32) \u{x)\ < Cd<?R l - n l 2 -*>\\Du\\ LHnRm , 

where C = C(n, m, v,jUq, Ni). 
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Proof. Let u be an extension of u by zero on Br(x) \ Q. Notice that u e W 1 ' 2 (Br(x)) and 
Du = xn R Du in Br(x). Then by the Poincare's inequality and (TLHb. we find that for all 
r E (0,/?/2] andy e Br/2(x), we have 

f |fi - fi r | 2 < Cr 2 \ \Du\ 2 = Cr \ \Du\ 2 

<C^R-" +2 - 2 ^\\Duf Ll{(lR{x)y 
Then by the Campanato's characterization of Holder seminorms, we have 

(5.33) [«]c*>ciM»» ^ CR x -" l2 -^\\Du\\ L 2 (aR(x)) . 

For any r e (d x ,R/2), there is x' e Br/2(x) \ Q. such that |x - x'| = r. By ( 15.331 ) we obtain 

|«(x)| = |«(x)-«(x')| < C^R l - nl2 -^\\Du\\ maR u)y 
By taking limit r — > <i A in the above inequality, we derive ( 15.321 i. ■ 

Now we are ready to prove the claim ( 15.30b . We may assume that d x < \x—y\/A because 
otherwise ( 15.301 ) follows from ( 15.291 ). We then set R = \x — y\/2 and u to be £-th column of 
G(-,y), for k — 1, . . . ,m, in Lemma [5.31| to obtain 

\G(x,y)\ < Cd^R l - nl2 -^\\DG{-,y)\\ L 2 (nK(x)) - R = \x- y\/2. 

On the other hand, since £Ir(x) cQ \ and R < R max /2, we have by ( 13.8b that 

\\DG(;y)\\ L 2 {nR{x)) < \\DG(;y)\y ( ^ BB(y)) < CR (2 - n)j2 . 

By combining the above two inequalities, we find that 

\G(x,y)\ < Cd'i \x -y\ 2 -"-'"', 

which implies ( 15.301 since we assume d x < \x - y|/4. We have proved the claim. 

We prove the estimate ( 13.14b using ( 15.301 ). Since the condition (1LHI) is symmetric be- 
tween L and 'L, by applying the above argument to 'G(x, y) and then interchanging x and y, 
we obtain, via the identity ( 13.7b and Remark [3.121 that 

(5.34) \G(x,y)\ < C{d y A |jc - y\f°\x - yt"-^ if < \x - y\ < 2R max . 

Again, we may assume that d x < \x — y|/8 to prove (13. 14b because otherwise ( 13.14b 
would follow from ( 15.341 ). We set R = \x — y|/4 and u to be A:-th column of G(-,y), for 
k — 1, . . . ,m, in Lemma l5.3U and then use the Caccioppoli's inequality to obtain 

(5.35) \G(x,y)\ < CdTR l -" l2 -^\\DG{-,y)\\ m a R(x)) < Cd^R-" /2 -^\\G(-,y)\\ L 2 (n2R{x)) . 

Notice that for all z e Q.2r{x), we have 2R < \z - y\ < 6R. Therefore, by the assumption 
R = \ x - y\/4 and ( 15.34b . we obtain 

(5.36) \G(z,y)\<C{d y A\x-y\r\x -y\ 2 -"-K\ Vz 6 Q 2R (x). 
By combining ( 15.351 ) and ( 15.361 ). we obtain 

\G(x,y)\ < Cd^x-yr^ldy A \x - y\f°\x - y^"^, 

which implies ( 13.141 ) since we assume d x < \x - y|/8. This completes the proof of ( 13.14b 
for all x,y E Q. satisfying < \x - y\ < R max . 

Next, we prove the second part of the theorem. Suppose R max < oo and diam(Q) < oo. 
Let x,y be arbitrary but fixed points in Q satisfying \x -y\ > R max /2. Let R = R ma xl^ and v 
be the k-Xh column of G(-,y) for k — 1, . . . , m. Notice that v e W 1 ' 2 (Q.r(x)) and v is a weak 
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solution of Lv =0 in Br(x) vanishing on Zr{x). Hence, by ( 15.61 with p = 2n/(n - 2), we 
have 

Therefore, by the above estimate and ( 13.8b we have 

(5.37) \G(x,y)\ < Cfl (2 -" )/2 ||G(-,y)|| y ,, (n ^ s(v)) < CR 2 " < CR 2 ma " x . 

On the other hand, if we set R = R max /4 and u to be the £-th column of G(-,y) for 
k — 1, . . . ,m, in Lemma IB. 311 then by ( 13.81 ) again, we find that if d x < R/4- = R max /8, then 

\G(x,y)\ < Cd^R^-^DGi-^h^)) < Cd^R 2 ^. 
By combining ( 15.371 ) and the above estimate, we derive the following conclusion. 

(5.38) \G(x, y)| < C( £ /, A R^TR 2 ^ whenever |x - y| > 
Then, by using (15.381 ) and arguing similarly as above, we obtain 

(5.39) \G(x, 3OI < C(d x A R mnx r°(d y A R max f° R 2 J; 2tla whenever |* - y| > fl„ M ,. 
Therefore, we conclude from ( 15.391 ) that for all x, y e Q satisfying \x - y\ > R max , we have 

\G(x,y)\ < C{d x A |x — A |jc - yir^max/ diamfi) 2 -"- 2 ^^ - y \ 2 -"- 2 f\ 

From the above estimate, we obtain (13.141 ) in case when \x — y\ > R max , with the constant C 
replaced by (R max l diam£T) 2 ~"~ 2w, C. Recall that we already have (13.14b in the case when 
< \x - y\ < Rmax- The proof is complete. ■ 

6. Appendix 

Lemma 6.1. Assume the condition (ILHb . For any p e (n/2, oa], there exists a constant 
C — C(n,?n,v,/j.Q,Ni,p) such that for all x e Q, R e (0,R max ), and f e L p (£Ir(x)), the 
following holds: If u € W'' 2 (Q s (x)) is a weak solution of either Lu — f or 'Lu — f in 
£Ir(x) and vanishes on Sr(x), then we have 

(6.2) Nli.»(n 8/2 ) < C (*~" /2 ||«|b ( n s) + R 2 -" ,l, \\f\\L^n R} ) ; O s = Cl R (x). 

In particular, the condition (ILBb holds with the same R max and No — No(n, m, v,po,Ni). 

Proof. We shall only consider the case when u is a weak solution of Lu — f since the 
proof of the other case is identical. Throughout the proof, we denote by C a constant 
depending on the prescribed parameters n, m, v, p and also the numbers fio,Ni that appears 
in the condition (1LHI ). As usual, the constant C may vary from line to line. 

Fix R < R max /4- and let u be a weak solution of Lu — f in Q4R = Q.4 R (xq) vanishing 
on E4R, where / e L p (D.4r) with p e (n/2, 00]. Fix x e Qr and s e (0,7?]. We write 
u — v + w in Q. s (x), where v € W l - 2 (Q. s (x)) is a weak solution of Lv — in Q. s (x) such that 
v -u € W ( J' 2 (f2. s (x)). Notice that v vanishes on E s {x). Then, (1LH1 ) implies that for < r < s, 

f |Dm| 2 < 2 f |Dv| 2 + 2 f |Dm>| 2 

Jn r (x) Jn,w 

<c(-)"" +2/, ° f \Dv\ 2 +2 f \Dw\ 2 
< c(-)" " +2/i ° f |D«| 2 + C f |£>h>| 2 . 
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Observe that w e W^' 2 (Q. s (x)) and w is a weak solution of Lw = f in Q. s (x). Therefore, we 
obtain 

f \Dw\ 2 < Q\f\^ 
Choose po e (n/2,p) such that jji := 2 - n/po < A'o- Then 

HflfU*>«. W) < ll/llL™ ( n lW )l". S | 1+2/ '" 2/P0 * C||/||i, 0(n2fl) ,"- 2+2 ^. 
By combining the above inequalities, we have for all r < s < R 



f \Du\ 2 < C (-)" 2+ ^° f \Du\ 2 + Cs"- 



A well known iteration argument (see e.g., (9J §111.2]) yields that for all r e [Q,R] and 
x e Q«, we have 

\Du\ 2 <C - \Du\ 2 + Cr"- 2 

Let u be an extension of u by zero on Z?2« \ Qir- Notice that u e W x ' 2 (B2r) and 
Du = xn^ R Du in Z?2«- Then by the Poincare's inequality and ( 16.3l l. we find that for all 
r e (0, /?] and x e Br, we have 

(6-4) f |fi - h,| 2 < Cr" +2 ^' («-" +2 -^' ||DH|| 2 2(nifi) + ||/|ll™ ( n 2 J ■ 

Then it follows from (16.4b and Holder's inequality that 

[«] 2 P1(Bs) < CR-"^\\Du\\l Hati!) + CJ^'-^ig,^. 

Therefore, we obtain 

ll«IIZ» ( n s/2) ^C^[«] 2 P1(Bfl) + C^"||«|| 2 2(Bs) 

< C^-" +2 ||D«|| 2 2(n2s) + C^-^II/II 2 ,^ + C/?-"||«|| 2 2(n;;) . 
Recall that u vanishes on Eur. By the Caccioppoli's inequality, we derive 

\\D< H n m) * CR- 2 \\uf L2(SliR) + C||/|| 2 2 „ /( „ t2)(n2fl) 

< c/?- 2 ii«n 2 2(n4s) + c^ 2+ "- 2 "^ii/ni„ (n4s) . 

By combining the above two inequalities and replacing R by R/4-, we obtain 

IMlL~(n s/8 ) < CR- n/2 \\u\\ LH n s) + CR 2 - n/ "\\f\\ m n R) . 

Finally, the above inequality together with a standard covering argument yields ( 16.2b . The 
proof is complete. ■ 

Acknowledgment. We thank the referee for useful comments. Kyungkeun Kang was 
supported by the Korean Research Foundation Grant (MOEHRD, Basic Research Promo- 
tion Fund, KRF-2008-33 1-C00024) and the National Research Foundation of Korea(NRF) 
funded by the Ministry of Education, Science and Technology (2009-0088692). Kyungkeun 
Kang appreciates the hospitality of the Department of Computational Science and Engi- 
neering, Yonsei University. Seick Kim was supported by the Korea Science and Engi- 
neering Foundation grant (MEST, No. R01-2008-000-20010-0) and also by WCU(World 
Class University) program through the Korea Science and Engineering Foundation (MEST, 
No. R31-2008-000-10049-0). 



GREEN'S MATRIX 



17 



References 

[I] Alfonseca, A.; Auscher, P.; Axelsson, A.; Hofmann, S.; Kim, S. Analyticity of layer potentials and L? 
Solvability of boundary value problems for divergence form elliptic equations with complex L°° coefficients. 
larXiv:0705.0836V l [math.AP] 

[2] Auscher, P.; Tchamitchian, Ph. Gaussian estimates for second order elliptic divergence operators on Lipschitz 

and C' domains. Evolution equations and their applications in physical and life sciences (Bad Herrenalb, 1998), 

15-32, Lecture Notes in Pure and Appl. Math., 215, Dekker, New York, 2001. 
[3] Chang, T.; Choe, H. J. Estimates of the Green 's functions for the elasto-static equations and Stokes equations 

in a three dimensional Lipschitz domain. Potential Anal. 30 (2009), no. 1, 85-99. 
[4] Cho, S.; Dong, H.; Kim, S. Global estimates for Green's matrix of second order parabolic systems with 

application to elliptic systems in two dimensional domains, preprint. 

[5] Dong, H.; Kim, S. Green's matrices of second order elliptic systems with measurable coefficients in two 
dimensional domains. Trans. Amer. Math. Soc. 361 (2009), no. 6, 3303-3323. 

[6] Dolzmann, G.; Muller, S. Estimates for Green's matrices of elliptic systems by L 1 ' theory. Manuscripta Math. 
88 (1995), no. 2, 261-273. 

[7] Fuchs, M. The Green-matrix for elliptic systems which satisfy the Legendre-Hadamard condition. 
Manuscripta Math. 46 (1984), no. 1-3, 97-115. 

[8] Fuchs, M. The Green matrix for strongly elliptic systems of second order with continuous coefficients. Z. 
Anal. Anwendungen 5 (1986), no. 6, 507-531. 

[9] Giaquinta, M. Multiple integrals in the calculus of variations and nonlinear elliptic systems. Princeton Uni- 
versity Press, Princeton, NJ, 1983. 

[10] Giaquinta, M. Introduction to regularity theory for nonlinear elliptic systems. Birkhauser Verlag, Basel, 
1993. 

[II] Gilbarg, D.; Trudinger, N. S. Elliptic partial differential equations of second order. Reprint of the 1998 ed. 
Springer- Verlag, Berlin, 2001. 

[12] Griiter, M.; Widman, K.-O. The Green function for uniformly elliptic equations. Manuscripta Math. 37 
(1982), no. 3,303-342. 

[13] Hofmann, S.; Kim, S. The Green function estimates for strongly elliptic systems of second order. 
Manuscripta Math. 124 (2007), no. 2, 139-172. 

[14] Littman, W.; Stampacchia, G; Weinberger, H. F. Regular points for elliptic equations with discontinuous 
coefficients. Ann. Scuola Norm. Sup. Pisa (3) 17 (1963), 43-77. 

[15] Maly, J.; Ziemer, W. P. Fine regularity of solutions of elliptic partial differential equations. American Math- 
ematical Society, Providence, RI, 1997. 

[16] Sarason, D. Functions of vanishing mean oscillation. Trans. Amer. Math. Soc. 207 (1975), 391^105. 
[17] Shen, Z. A note on the Dirichlet problem for the Stokes system in Lipschitz domains. Proc. Amer. Math. Soc. 
123(1995), no. 3, 801-811. 

(K. Kang) Department of Mathematics, Sungkyunkwan University, Suwon 440-746, Republic of Korea 
E-mail address: kkang@skku . edu 

(S. Kim) Department of Mathematics, Yonsei University, Seoul 120-749, Republic of Korea 
Current address: Department of Computational Science and Engineering, Yonsei University, Seoul 120-749, 
Republic of Korea 

E-mail address: kimseickOyonsei . ac . kr 



